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电子 断层 三 维 重 构 技术 的 研究 与 进展 


张 法 ， 万 晓 华 ， 储 琪 


摘要 : 电子 断层 三 维 重 构 技 术 能 够 在 纳米 尺度 下 重 构 出 不 具有 全 同性 的 细胞 或 大 分 子 的 三 维 结构 ， 已 经 成 
为 一 种 公认 的 研究 生物 大 分 子 结构 强 有 力 手 段 。 然 而 ， 电 子 断 层 三 维 重 构 仍 然 存 在 二 维 图 像 对 位 不 精确 、 
重 构 算法 精度 低 和 重 构 速度 极其 绥 慢 等 问题 。 针 对 以 上 问题 ， 本 文 详细 介绍 了 电子 断层 三 维 重 构 的 发 展 和 
现状 ， 分 析 了 当前 所 面临 的 主要 挑战 问题 ， 着 重 介绍 了 我 们 在 电子 断层 三 维 重 构 相 关 研 究 工 作 上 的 进展 。 


起 


关键 词 : 电子 断层 三 维 重 构 对 位 算法 迭代 算法 并 行 计 算 GPU 


~ 结构 生物 学 是 研究 生物 大 分 子 特定 三 维 结构 及 结构 的 特定 运动 、 结 构 相 关 生 物 学 功能 的 
T 有 一 门 学 科 ， 其 起 源 可 追溯 到 上 世纪 50 年 代 沃 和 森 (J. Waston) 和 克 里 克 〈F. Crick) 发 现 DNA 
> 双 螺 旋 结构 ,60 年 代 佩 鲁 茨 (M. Perutz) 利 用 X 射 线 唱 体 衍 射 技术 获得 肌 球 蛋白 的 三 维 结构 口 。 
三 现在 结构 生物 学 在 分 子 生物 学 研究 中 占据 了 主流 地 位 , 已 成 为 药物 设计 、 疫苗 开发 等 研究 的 
基础 。 


众所周知 ,蛋白质 是 生命 活动 的 体现 者 之 一 , 对 重 白 质 三 维 结构 的 研究 一 直 是 结构 生物 

学 的 主要 研究 内 容 。 目 前 X 射线 晶体 衍射 技术 、 核 磁 共 振 技术 和 电镜 三 维 重 构 技术 是 结构 

生物 学 的 三 大 研究 手段 。 其 中 X 射线 晶体 衍射 技术 是 目前 分 辨 率 最 高 的 结构 测定 方法 ， 可 

以 获得 生物 分 子 的 原子 分 辩 率 三 维 结构 , 但 是 生物 大 分 子 的 晶体 难以 获得 , 结构 解析 极其 困 

> 难 ， 核 磁 共 振 技术 可 以 获得 和 蛋白质 在 溶液 中 的 三 维 结构 ， 能 够 获得 蛋白 质 的 动态 变化 信息 ， 

SZ 但 是 其 研究 对 象 的 分 子 量 通 常 难以 超过 20kD， 因 此 这 两 种 方法 ， 极 大 限制 了 人 类 对 蛋白 质 

©& 功能 以 及 生命 活动 规律 的 认识 。 随 着 生物 样品 制备 技术 的 不 断 完善 , 电子 显微镜 设备 的 进步 

~ 以 及 数字 图 像 处 理 技术 的 发 展 ， 电 镜 三 维 重 构 技 术 ， 特 别 是 冷冻 电镜 三 维 重 构 技 术 ， 能 够 解 

析 那 些 无 法 应 用 X 射线 晶体 衍射 和 核磁 共振 技术 进行 分 析 的 蛋白 质 ， 弥 补 这 两 种 技术 的 不 

是 ， 已 成 为 研究 生物 大 分 子 三 维 结构 的 重要 手段 。 萨 利 (A. Sali) 等 在 《Nature》 提 出 了 今 

后 结构 生物 学 研究 的 发 展 方向 ， 即 将 包括 X 射线 晶体 衍 财 、 核 磁 共 振 、 冷 冻 电 镜 三 维 重 构 、 

菊 光 能 量 共振 等 多 种 生物 物理 技术 联合 起 来 , 在 分 子 水 平 上 理解 细胞 的 结构 ,其 中 冷冻 电镜 
三 维 重 构 技术 起 到 了 重要 的 纽带 作用 中 。 


E 镜 三 维 重 构 技 术 的 思想 最 早 由 德 罗 西 耶 (D. DeRosier) 和 克 卢 格 (A. Klug) 在 1968 
年 提出 。 上 世纪 80 年 代 ， 快 速 冷冻 生物 样品 的 制备 技术 和 低温 电子 显微镜 技术 的 快速 发 展 
与 实用 化 ， 为 利用 电镜 三 维 重 构 技 术 研究 和 蛋白质 三 维 结构 黄 定 了 基础 。 当 前 ， 利 用 电镜 三 维 
重 构 技 术 解 析 蛋 白质 超 分 子 复 合 物 的 结构 已 分 化 发 展 为 三 种 具有 不 同 特点 和 适用 范围 的 方 
法 : 电子 唱 体 学 (Electron Crystallography)、 单 颗粒 三 维 重 构 (Single Particle) 和 电子 断层 
三 维 重 构 (Electron Tomography，ET)。 其 中 ， 电 子 唱 体 学 主要 用 于 解析 和 蛋白质 分 子 的 二 维 
晶体 样品 ;， 单 颗粒 三 维 重 构 ， 主 要 解析 具有 全 同性 的 蛋白 质 分 子 的 三 维 结构 ， 最 新 研究 成 果 
已 经 可 以 在 原子 分 辩 率 下 解析 具有 高 对 称 性 的 病毒 分 子 三 维 结构 (关于 单 颗粒 三 维 重 构 的 
详细 介绍 ， 参 见 本 刊 第 56 期 “ 单 颗粒 冷冻 电镜 三 维 重 构 研 究 的 新 进展 ”一 文 由 )， 电 子 断层 
三 维 重 构 通 过 获取 同一 区 域 多 个 角度 的 投影 图 来 反问 重 构 所 研究 对 象 的 三 维 结构 , 适合 于 在 
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获得 结构 的 分 辨 率 〈 约 4~10 纳米 ) 不 能 与 
同性 生物 样品 的 三 维 结构 和 功能 


定形 、 不 对 称 和 不 具有 全 


电子 断层 三 维 


结构 。 由 于 无 需 分 子 样品 具有 结构 同一 性 和 对 称 性 ， 因 此， 虽然 目前 电 


电子 断层 三 维 重 构 的 适用 


斥 度 非常 广泛 ,包括 从 分 子 水 平 的 蛋 
以 至 细胞 水 平 的 组 织 结 构 ， 可 以 有 效 弥 补 X 射线 晶体 衍射 、 核 磁 
方法 得 到 的 高 精度 结构 与 光学 显 
解 生物 大 分 子 的 功能 提供 重要 而 有 益 的 信息 。 因 此 ， 电 子 断 层 三 允 


微 技术 得 到 的 低 分 辩 率 纪 


构 技 术 的 研究 与 进 必 


纳米 级 尺度 上 研究 不 具有 结构 均一 性 的 蛋白 、 病毒 、 细 胞 器 以 及 它们 之 间 组 成 的 复合 体 的 三 
多 


子 断 层 三 维 重 构 所 


有 着 不 可 替代 


其 它 两 种 电镜 三 维 重 构 方法 相 比 ， 但 其 在 研究 非 


的 重要 作用 。 此 外 ， 


白质 , 到 亚 细 胞 水 平 的 细胞 器 ， 


* 振 、 单 颗粒 三 维 重 构 等 
j 胞 整体 图 像 之 间 的 空白 , 为 深入 理 
重 构 技术 被 《科学 》 和 杂志 


评 为 2002 年 十 大 科技 突破 之 一 四， 近年 来 在 冷冻 电镜 研究 领域 得 到 了 越 来 越 多 的 关注 。 


2 ”电子 断层 三 维 重 构 技 术 


电子 断层 三 维 重 构 技术 (ET) 是 从 一 个 物体 的 电子 显 微 投影 图 像 重 构 物体 内 部 结构 的 
技术 。 类 似 于 人 们 熟知 的 CT 扫描 技术 ， 通 过 获取 同一 个 物体 的 多 个 连续 角度 下 的 二 维 投影 


图 像 来 反 向 重 构 它 的 三 维 结构 。 其 主要 策 


让 样品 沿 着 


个 与 电子 束 垂直 的 轴 


二 维 投影 图 像 ， 然 后 对 这 一 系列 二 双 
结构 。 和 CT 扫描 技术 转动 光源 获得 不 同 / 
通过 倾 转生 物 样品 来 得 到 不 同方 向 的 投影 。 


TS 


+60°~+80° 之 间 ， 间隔 角度 为 1 或 2"。 因 此 ，ET 数据 存在 信息 缺失 ， 即 缺 


wedge)。 电 子 断 层 三 维 重 构 过 程 主要 包括 以 下 几 个 步骤 中 ，(D, 样品 


各 是 将 制备 好 的 生物 样品 放 在 透射 电子 显微镜 中 ， 
旋转 ,每 旋转 一 定 角度 ,拍摄 一 张 样品 在 相对 应 方向 上 的 
投影 图 像 进行 对 位 , 利用 三 维 重 构 算法 得 到 样品 的 三 维 
和 度 的 投影 图 不 同 ，ET 的 光源 〈 电 子 束 ) 固定 ， 
根据 电镜 设备 的 不 同 ， 样 品 倾 转 范围 大 概 在 


失 栋 问题 (missing 


| 备 、(2). 获取 二 维 投 


影 图 像 、(3). 投影 图 像 的 对 位 、(4). 断层 图 像 的 三 维 重 构 、(5). 重 构 结果 的 后 期 处 理 ， 如 图 


1 所 示 叫 。 


样品 制备 : 电子 断层 重 构 技 术 主 要 


研究 细胞 环境 下 的 三 维 结构 , 与 传统 的 
CT 扫描 使 用 的 X 射 线 相 比 , 电子 束 能 罕 


透 的 样品 厚度 较 小 , 限制 


了 电子 断层 技 


术 对 大 的 生物 样品 (如 完整 的 真 核 细 
胞 ) 的 研究 。 对 于 电子 束 难 以 穿 透 的 厚 


样品 ,如 细胞 等 ， 需 要 采用 切片 的 方法 


进行 制 样 ， 经 过 脱水 、 染 


色 后 ， 进 行 观 


察 和 拍摄 对 于 电子 束 能 
较 薄 的 样品 ， 如 病毒 、 细 
用 快速 冷冻 的 方法 准备 术 


穿 透 的 较 小 或 
菌 等 ,可 以 采 
lf 品 并 进行 冷 


冻 电 子 断 层 三 维 重 构 。 冷 冻 电 子 断 层 三 


状态 , 克服 了 由 于 脱水 和 
来 的 假象 , 较 好 地 保护 了 
的 完整 性 。 


维 重 构 的 主要 优势 是 样品 接近 于 天 然 


染色 等 可 能 带 
生物 样品 结构 


二 维 投影 图 像 的 获取 : 电子 断层 的 
二 维 投影 图 像 收集 首先 是 对 样品 选取 合适 的 成 像 区 域 , 然后 对 该 区 域 连续 收集 不 同 角度 的 成 


像 投影 照片 。 电 子 断 层 二 维 投影 区 
样品 成 像 区 域 的 电子 辐射 损伤 ， 一 般 利 月 


电子 检 
wl 原始 图 像 


图 1. 


ci 2 
4 


充 
重 构 
电子 断层 三 维 重 构 流程 


已 对 准 图 像 


= i 
示意 图 。 
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像 的 收集 过 程 大 概 如 下 : (1). 寻找 合适 区 域 , 为 了 减少 对 


日 低 电子 束 剂量 模式 的 数据 采集 方案 ; (2). 调节 样品 


第 9 卷 第 5 期 言 恩 技术 快报 Vol.9 No.5 

Information Technology Letter Sep. 2011 
台 旋 转 中 心 高 度 ， 使 样品 台 旋 转 时 样品 维持 在 成 像 中 心 范围 ; (3). 调 焦 ，(4). 找到 同一 成 像 
区 域 的 精确 跟踪 ; (5). 照相 ; (6). 样品 倾斜 到 下 一 个 角度 ， 重 新 找到 样品 的 成 像 区 域 中 心 ， 
然后 重复 以 上 四 步 : 寻找 、 调 焦 、 跟 踪 和 照相 。 


由 于 电子 断层 在 收集 二 维 投影 图 像 时 ， 需 要 对 样品 同一 区 域 反 复 成 像 , 往往 容易 造成 较 
大 幅度 的 辐 照 损伤 ， 因 此 必须 严格 控制 样品 台 的 精确 倾 转 ， 对 成 像 区 域 的 自动 跟踪 ， 以 及 每 
次 成 像 的 电子 束 剂 量 等 因素 。 在 获取 投影 图 像 的 过 程 中 ,样品 在 旋转 时 由 于 机 械 原 因 ， 往 往 
会 有 一 定 的 偏 移 , 需要 对 不 同 倾斜 角度 的 图 像 通 过 互相 关 运 算 来 完成 精确 跟踪 ,为 了 提高 图 
像 衬 度 和 跟踪 的 准确 性 ,可 以 在 样品 中 加 入 胶体 金 颗粒 。 目 前 已 有 一 些 自动 化 的 电子 断层 数 
据 收集 软件 ， 比 如 UCSF tomographyl"，FEI 公司 的 Explor3D 等 。 这 些 软件 能 够 利用 最 初 几 
张 图 像 的 偏 移 算出 整个 数据 收集 过 程 中 图 像 的 储 移 轨迹 , 从 而 大 大 降低 了 额外 增加 的 剂量 对 
样品 的 损伤 ， 并 且 能 够 比较 精确 地 控制 电镜 和 CCD 相机 ， 完 成 自动 聚集 功能 。 


二 维 投影 图 像 的 对 位 : 由 于 在 收集 二 维 投影 图 像 时 ， 样 品 在 倾斜 过 程 中 不 可 避免 地 会 出 
现 倾 转 轴 改 变 、 样 品 本 身 移 动 、 投 影 本 身 移 动 、 投 影 拍摄 时 CCD 照 相机 参数 不 一 致 等 问题 ， 
使 得 投影 的 系列 二 维 图 像 之 间 发 生 偏 移 、 旋 转 或 者 放大 缩小 等 效应 。 因 此 , 在 三 维 重 构 之 前 ， 
~ 必须 对 一 个 倾斜 系列 投影 的 图 像 进 行 对 位 , 首先 是 计算 倾斜 过 程 中 由 于 样品 台 等 的 移动 造成 
F 的 图 像 中 心 的 移动 ; 其 次 , 还 需要 对 由 于 现存 的 成 像 条 件 导 致 的 倾斜 角度 , 图 像 本 身 的 旋转 、 
扭曲 以 及 放大 倍数 等 进行 修正 。 二 维 投影 图 像 对 位 的 质量 直接 决定 了 最 终 重 构 结果 的 质量 。 


目前 , 二 维 投 影 图 像 对 位 常用 的 方法 主要 有 基于 标记 点 的 方法 和 基于 图 像 互 相关 匹配 的 
方法 。 基 于 标记 点 的 方法 是 在 样品 制备 过 程 中 , 加 入 若干 一 定 大 小 的 胶体 金 (golden particles ) 
作为 标记 物 ， 然 后 根据 胶体 金 标记 的 对 应 ， 建 立 投 影 模型 ， 求 解 偏 移 、 旋 转 或 放大 缩小 等 对 
位 参数 。 这 种 方法 常用 的 软件 是 IMOD 色 。 但 这 种 方法 所 添加 的 胶体 金 极 有 可 能 造成 样品 污 
染 ， 而 且 在 倾 转 过 程 中 胶体 金 可 能 会 发 生 漂移 ， 导 致 无 法 很 好 地 对 位 。 另 一 种 对 位 方法 是 利 
用 图 像 本 身 的 互相 关 分 析 对 图 像 进行 对 位 ， 其 中 比较 常用 的 软件 有 Protomog。 除 了 IMOD 
和 Protomo 外 , 还 有 其 他 一 些 软件 用 于 二 维 图 像 的 对 位 ， 如 Spider 等 以 及 FEI 的 Inspec3D 
等 商业 软件 ， 但 是 精确 的 二 维 图 像 对 位 一 直 是 该 领域 的 一 个 难点 问题 。 


电子 断层 三 维 重 构 : 电子 断层 的 三 维 重 构 是 由 所 有 对 位 的 二 维 投影 图 像 重 新 反 向 投影 到 
一 个 3D 体 积 中 形成 的 ， 其 算法 主要 有 两 种 : 背 投影 重 构 和 迭代 重 构 。 


9 背 投影 重 构 算 法 包括 直接 背 投影 (back projection, BP)、 滤 波 背 投影 (filter back projection, 
FBP) 和 加 权 背 投影 (weighted back projection，WBP)。 直 接 背 投影 算法 的 基本 思想 是 三 维 
物体 某 一 点 的 密度 等 于 通过 该 点 所 有 射线 密度 的 总 和 。 直 接 背 投影 算法 是 一 种 最 简单 的 重 构 
算法 , 但 它 容易 使 重 构 结果 在 实 空间 产生 点 扩散 的 不 良 效果 。 滤波 背 投 影 算法 提出 对 重 构 结 
果 进 行 点 扩散 函数 的 去 卷 积 处 理 , 但 由 于 样品 的 空间 几何 参数 未 知 , 无 法 得 到 有 效 的 点 扩散 
函数 。 加 权 背 投影 算法 则 是 对 二 维 投影 图 像 传 里 时 变换 的 不 同 频 率 进 行 加 权 处 理 ， 然 后 在 实 
空间 中 进行 反 投影 操作 。 迭代 重 构 算 法 完全 是 在 实 空间 进行 ,建立 三 维 物体 中 所 有 点 的 密度 
与 投影 图 中 所 有 点 的 密度 之 间 的 线性 方程 组 ， 通 过 迭代 方法 求解 线性 方程 组 得 到 三 维 模型 ， 
使 得 原始 投影 图 和 重 构 投影 图 尽量 接近 。 常 用 的 迭代 重 构 算法 有 代数 重 构 算法 (Algebraic 
Reconstruction Technique, ART)、 同步 代 数 重 构 算 法 上 (Simultaneous Algebraic Reconstruction 
Technique，SART) 和 联合 迭代 重 构 算 法 0 (Simultaneous Iterative reconstruction Technique， 
SIRT)。 友 代 重 构 算法 具有 一 致 的 权 函 数 ， 但 计算 量 巨大 ， 为 了 得 到 满意 的 重 构 结果 往往 需 
要 花费 大 量 时 间 。 


重 构 结 果 的 后 期 处 理 : 为 了 减少 对 样品 的 损伤 ， 对 样品 拍照 时 往往 电子 束 剂 量 很 少 ， 所 


一 轴 
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以 得 到 的 二 维 投影 图 像 的 信 噪 比 都 比较 低 ， 由 此 产生 的 重 构 结构 噪 音 也 很 大 。 所 以 对 重建 结 
果 要 进行 适当 的 降 噪 处理 以 突出 研究 对 象 。 常 用 的 去 噪 方法 包括 中 值 滤波 、 非 线性 滤波 等 。 
其 中 , 各 向 异性 扩散 的 非 线性 滤波 方法 效果 较 好 。 为 了 进一步 将 研究 对 象 从 背景 和 噪声 中 区 
分 开 ， 需 要 对 重建 结果 进行 分 割 。 最 常用 的 方法 是 使 用 手动 分 制 ， 需 要 用 户 绘制 对 象 和 背景 
的 边界 轮廓 。 分 割 完成 后 可 以 利用 表面 泻 染 技术 对 重 构 结果 进行 可 视 化 。 常 用 的 分 割 软件 包 
括 IMOD 四 和 VAT4ML31。 


2.1 当前 研究 面临 的 问题 


近年 来 电子 断层 技术 取得 了 长 足 的 进展 , 已 经 成 为 结构 生物 学 研究 中 至 关 重 要 的 研究 手 
段 。 我 国 在 电子 断层 重 构 领域 的 研究 也 取得 了 迅猛 的 进展 。 国 内 已 拥有 了 包括 2 台 Titan 
Krios、1 台 Tecnai F30 等 在 内 的 高 端 电子 显微镜 研究 设备 。 同 时 ， 多 家 研究 机 构 也 计划 或 正 
在 购置 更 多 的 电镜 设备 。 国 内 已 经 具备 了 很 好 的 契机 来 开展 电子 断层 成 像 方面 的 工作 。 然 而 
与 迅猛 发 展 的 硬件 设施 相 比 , 我 们 在 相应 软件 方面 的 研究 就 显得 严重 滞后 , 存在 以 下 三 方面 
的 问题 ， 制 约 着 电子 断层 重 构 技术 的 进一步 发 展 : 


EE 1. 缺乏 准确 的 二 维 投影 图 像 对 位 算法 


如 前 所 述 , 在 电子 断层 成 像 对 样品 投影 成 像 过 程 中 , 各 种 因素 会 导致 最 终生 成 的 图 像 序 
列 之 间 出 现 偏 移 、 旋 转 或 者 放大 缩小 等 效应 , 所 以 在 三 维 重 构 前 需要 对 系列 二 维 投影 图 像 进 
行 对 位 。 

目前 , 最 常用 的 对 位 方法 是 前 述 之 胶体 金 标记 法 和 无 标记 对 位 方法 。 胶体 金 标 记 法 存在 
很 多 缺点 ,所 以 需要 更 普 适 的 无 标记 方法 来 对 位 投影 图 像 。 无 标记 对 位 方法 主要 有 以 下 几 类 : 
(D. 互相 关 方法 , 其 缺陷 在 于 不 同 角 度 下 的 投影 成 像 有 可 能 差别 很 大 , 计算 的 偏 移 值 不 准确 ， 
且 容 易 造成 累积 误差 (2). 等 价 线 方法 ， 其 局 限 性 在 于 上 只 能 求解 一 个 对 位 参数 ， 无 法 准确 地 
对 位 投影 图 像 ，(3). 迭代 对 位 方法 ， 需 要 进行 一 次 三 维 重 构 和 求解 参数 进行 对 位 操作 ， 计 算 
量 巨 大 ; (4). 基于 特征 匹配 的 对 位 方法 ， 使 用 互相 关 方 法 来 匹配 特征 点 ， 存 在 着 由 于 互相 关 
方法 所 带 来 的 点 与 点 之 间 对 应 不 精确 的 问题 。 因 此 ， 电 子 断 层 技术 人 迫切 需要 更 为 精确 的 二 维 
投影 图 像 对 位 方法 。 


2. 缺乏 高 精度 的 三 维 重 构 算法 


电子 断层 三 维 重 构 算 法 主要 包括 加 权 背 投影 和 迭代 算法 两 种 。 加 权 背 投影 是 由 雷 德 麦 彻 
CM. Rademacher) 提出 的 最 早 的 电子 断层 重 构 算 法 。 它 利用 中 央 截 面 定理 ， 其 基本 原理 是 
将 每 个 二 维 投影 图 像 背 向 其 倾斜 角 方 向 投影 到 三 维 空间 中 , 多 个 背 投 影 图 像 在 三 维 空间 中 县 
加 形成 样品 的 三 维 结构 。 在 重 构 过 程 中 ,通过 权重 滤波 器 实现 取样 平均 的 目的 。 由 于 加 权 背 
投影 算法 简单 便于 实现 , 目前 流行 的 大 部 分 主流 电子 断层 重 构 软件 , 如 IMOD™, SPIDERO 
等 均 采用 加 权 背 投影 算法 进行 三 维 重 构 。 然 而 ， 加 权 背 投影 算法 容易 受到 信息 缺失 的 影响 ， 
导致 z 轴 方 向 上 的 重 构 结果 错误 , 此 外 加 权 背 投影 算法 在 重 构 过 程 中 没有 考虑 噪声 带 来 的 影 
响 ， 使 得 重 构 结果 不 精确 。 


近 几 年 人 们 开始 考虑 采用 从 代 算 法 进行 电子 断层 三 维 重 构 。 由 于 电子 断层 中 投影 方向 只 
能 覆盖 有 限 角 范 围 , 并 且 投 影 图片 噪 声 很 大 , 迭代 算法 在 这 种 数据 不 足 和 高 噪声 的 条 件 下 能 
够 获得 较 好 的 重 构 结 果 。 然 而 ， 传 统 的 迭代 算法 具有 以 下 缺陷 : (1). 往往 任意 选取 重 构 像 素 
值 的 初 值 ， 通 常情 况 下 选取 全 0 《或 全 1) 向 量 ， 可 能 偏离 真实 值 较 远 ， 因 此 需要 多 次 达 代 
才能 达到 所 需 的 收敛 值 ， (2). 在 更 新 投影 数据 时 采用 的 访问 方式 是 顺序 访问 策略 ， 由 于 相 
邻 投 影 方 向 之 间 有 很 大 的 相关 性 ， 容 易 引 入 误差 .影响 重 构 结果 的 精度 ， 并 且 降 低 达 代 算法 
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的 收敛 速度 ;” (3). 迷 代 算法 在 每 次 迭代 更 新 时 ， 误 差 的 补偿 只 考虑 加 权 和 矩阵 ， 这 显然 是 不 


合理 的 ， 因 为 投影 误差 是 由 加 权 和 矩阵 和 当前 像素 值 共 同 决定 的 。 
3. ”缺乏 高 效 的 三 维 重 构 并 行 算法 


随 着 电子 显 微 技 术 的 飞速 发 展 , 为 了 能 够 得 到 更 高 分 辨 率 的 重 构图 像 , 往往 需要 增加 投 
影 图 像 的 数量 和 尺寸 , 这 样 会 导致 电子 断层 三 维 重 构 需 要 大 量 的 计算 资源 和 处 理 时 间 。 然 而 ， 
面 对 这 些 大 数据 量 的 图 像 , 现 有 的 方法 无 法 进行 及 时 有 效 的 处 理 ， 已 经 成 为 限制 电子 断层 三 
维 重 构 广 泛 应 用 的 一 个 瓶颈 问题 。 围 绕 高 性 能 计算 技术 降低 三 维 重 构 时 间 ， 人 们 进行 了 大 量 
的 研究 2。 当前 大 部 分 电子 断层 三 维 重 构 软 件 ， 如 IMOD 包 ，SPIDERH0， 都 实现 了 传统 
并 行 计 算 架 构 〈 如 集群 环境 ) 上 的 算法 。 然 而 传统 集群 系统 的 昂贵 造价 以 及 GPU 硬件 技术 
的 迅猛 发 展 使 得 利用 GPU 开发 高 效 的 电子 断层 三 维 重 构 算 法 成 为 当前 研究 的 热点 02181。 


3 ”我 们 的 工作 


针对 以 上 述 目 前 电子 断层 存在 的 问题 , 本 节 主 要 介绍 我 们 课题 组 在 电子 断层 重 构 方面 的 
~ 一 些 研究 工作 。 首 先 ， 介 绍 在 电子 断层 中 一 种 新 的 基于 特征 提取 的 二 维 投影 图 像 对 位 算法 
F SAMA; 其 次 ,介绍 一 种 收敛 速度 更 快 ， 重 构 结 果 更 精确 的 迭代 重 构 算法 ASART; 最 后 ， 
介绍 了 我 们 开发 的 国内 首 款 基 于 GPU 平台 的 电子 断层 重 构 软 件 ATOM。 


3.1 二 维 投影 图 像 对 位 算法 SAMA 


针对 当前 无 胶体 金 样品 的 二 维 投影 图 像 对 位 问题 , 我 们 提出 了 一 种 基于 特征 匹配 的 无 标 
记 对 位 新 方法 一 SAMA(SIFT-based Automatic Markerless Alignment) 。 该 方法 建立 了 样品 在 投 
影 时 三 维 密度 点 与 投影 点 的 关系 的 非 线性 模型 , 于 是 对 位 问题 转化 为 根据 测量 样本 估计 模型 
参数 的 参数 估计 问题 , 通过 在 每 幅 二 维 投影 图 像 中 提取 特征 点 作为 测量 样本 点 , 根据 测量 样 
本 的 精确 性 来 选择 合适 的 非 线 性 模型 ， 并 进行 参数 估计 。 此 方法 的 框架 包括 : 特征 点 提取 、 
特征 点 匹配 、 特 征 点 跟踪 以 及 模型 选择 与 参数 求解 4 个 部 分 -SAMA 方法 的 具体 步骤 如 下 0 


特征 点 提取 : 图 像 的 特征 点 通常 包括 两 个 部 分 ，(1). 表示 子 (descriptor) 即 特征 点 独 
特 的 位 置 ， 一 般 根据 图 像 数 学 特征 进行 提取 ， 然 后 以 坐标 方式 显 式 表达 ; (2). 描述 子 


《description) 表示 子 领 域内 的 信息 ， 用 于 匹配 不 同 的 表示 子 。 早 前 的 特征 提取 方法 主要 关 
™ 注 于 如 何 设计 具有 区 分 度 的 表示 子 , 较 少 涉及 如 何 设计 具有 区 分 度 的 描述 子 。 我们 的 特征 点 


提取 采用 简化 的 尺度 不 变 特征 变换 (Scale Invariant Feature Transform，SIFT) 算 法 ， 对 表述 子 
和 描述 子 都 进行 加 强 。 其 中 表示 子 是 同一 图 像 不 同 尺 度 中 的 梯度 稳定 点 ; 描述 子 是 一 个 128 
维 向 量 来 描述 表示 子 周 围 灰 度 的 统计 信息 和 灰 度 方向 信息 。SIFT 算法 是 目前 最 为 精确 的 图 
像 特征 点 提取 方法 ， 简 化 的 SIFT 方法 能 显著 提高 其 速度 ， 但 不 影响 其 精确 性 。 


特征 点 匹配 : 由 于 图 像 上 的 每 个 特征 点 的 信息 都 可 以 用 一 个 128 维 向 量 描述 子 表示 , 且 
不 同 特征 点 的 描述 子 区 别 很 大 , 两 幅 图 像 之 间 可 通过 计算 最 小 欧式 距离 的 方法 进行 特征 点 匹 
配 。 此 外 两 幅 图 像 可 以 根据 相隔 角度 估计 出 一 幅 图 像 上 的 特征 点 在 男 一 幅 图 像 上 的 位 置 差别 
的 最 大 值 , 在 此 最 大 位 置 差别 区 域 之 内 搜索 与 之 欧式 距离 的 最 小 点 作为 匹配 特征 点 。 若 此 最 
小 欧式 距离 小 于 一 定 闵 值 , 则 将 此 匹配 记录 下 来 。 最 后 对 所 有 图 像 两 两 之 间 计 算 特 征 点 匹配 ， 
并 利用 两 幅 投 射 图 像 之 间 存 在 着 二 阶 几何 张 量 关 系 ， 去 除 错误 匹配 对 。 

特征 点 跟踪 : 通过 对 图 像 中 正确 的 匹配 对 进行 串联 来 形成 特征 链 。 每 条 特征 链表 示 同 一 
三 维 密度 点 的 投影 序列 。 初 始 时 , 对 一 幅 图 像 中 的 特征 点 a 找到 其 在 其 它 图 像 中 的 匹配 点 b， 
在 其 它 图 像 中 如 果 b 又 与 男 外 图 像 有 匹配 点 ce， 则 将 此 串联 起 来 。 这 样 对 每 个 图 像 对 ， 都 搜 
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索 原 有 的 特征 链 ， 若 匹配 对 对 应 在 原 有 某 条 特征 链 上 ， 则 此 链 上 加 点 ， 和 否则 重新 生成 新 特征 
链 。 


模型 选择 与 参数 估计 :求解 电子 断层 三 维 重 构 中 三 维 密度 点 与 二 维 投影 点 的 关系 模型 可 
表示 为 : 珊 ) = Ar +tj = sjR.(Q;)PR,(9,)R.(9;)r+t;。 我 们 需要 估计 s,Qj,9,,9jsrst; 这 6 个 
参数 。 其 中 声 ; 表示 第 j 幅 图 像 中 对 应 的 二 维 投影 点 坐标 ， 4 表示 从 三 维 坐 标 到 二 维 坐 标 
的 投影 矩阵 , rx; 表示 第 i 个 三 维 密度 点 , 5 为 成 像 时 二 维 投 影 点 的 平移 向 量 ，@ ,9,,9; 分 别 为 
绕 z，y,， x 轴 的 倾 转角 度 。 


在 参数 估计 过 程 中 我 们 使 用 捆绑 调整 (bundle adjustment) 和 重 棒 的 异常 点 (outlier〉 去 除 
方法 ， 根 据 误差 特性 对 参数 进行 估计 。 在 得 到 上 述 6 个 参数 后 ， 对 每 幅 图 像 进行 对 位 。 
使 用 SAMA 方法 对 真实 电子 断层 成 像 数 据 进行 对 位 , 参数 估计 残 差 值 在 0.2~0.4 像素 之 
间 ， 与 使 用 胶体 金 标 记 的 对 位 方法 进行 对 位 的 残 差 值 结果 “(0.2~0.7 像素 ) 具有 可 比 性 ， 且 
经 过 此 方法 对 位 后 重 构 的 图 像 与 使 用 胶体 金 方法 对 位 重 构 的 图 像 效 果 相当 。 

3.2 三 维 迭 代 重 构 算法 ASART 


迭代 算法 认为 重 构 的 三 维 结果 太 可 以 看 成 一 组 基 函 数 训 和 像素 灰 度 值 x 的 线性 组 合 , 表 
示 如 下 : 


J 
Pe ne 


其 中 , J 为 三 维 空间 像素 点 的 个 数 。 传 统 的 体 素 模 型 法 (voxel model) 以 单个 小 立方 体 ， 
即 体 素 (voxel)， 作 为 三 维 图 像 的 基本 单元 ， 使 用 的 基 函 数 为 分 段 函 数 ， 无 法 正确 地 体现 重 
构 样 品 密度 空间 的 连续 性 。 我 们 采用 一 种 blob 模型 ， 认 为 整个 三 维 空间 由 多 个 全 同 的 球体 
构成 (如 图 2 所 示 )， 由 于 球体 之 间 相 互 交 叉 重 登 ， 能 够 更 好 地 体现 密度 空间 的 连续 性 ， 使 
重 构 结果 更 加 精确 。 


根据 电子 断层 的 投影 原理 , 每 条 射线 的 投影 值 可 以 看 成 该 射线 经 过 的 像素 点 灰 度 值 的 加 
权 和 ， 用 如 下 线性 系统 表示 : 


> 
Pp; = WX ，， 
i ja 


其 中 ，p; 表 示 第 i 条 射线 的 投影 值 ，wj 表示 加 权 因 子 ， 
反映 了 像素 总 对 射线 投影 产 的 贡献 。 在 这 一 线性 模型 下 ， 电 
子 断 层 重 构 的 目的 是 采用 人 迭代 算法 根据 测量 所 得 的 产值 和 可 
以 计算 得 到 的 wj; 值 估计 各 个 像素 点 灰 度 值 x, 最 后 得 到 整个 
三 维 结构 f。 在 blob 模型 下 ， 由 于 每 个 球体 的 投影 值 跟 射线 
方向 无 关 ， 我 们 利用 一 种 印迹 (footprint) 方法 求解 wy， 先 
提前 计算 单个 球体 在 不 同 r 位 置 上 的 投影 值 并 存在 列表 ， 
然后 根据 像素 点 x 到 射线 pi 的 距离 查 表 得 到 对 应 的 wy 值 。 


在 求解 入 值 时 利用 现 有 的 和 迭代 算法 仍然 存在 收敛 速度 较 
慢 ， 结 果 不 精确 等 不 足 ， 因 此 我 们 提出 一 种 自 适应 联合 代数 
过 代 算法 ASART 0 ， 采 用 以 下 四 种 关键 技术 来 提高 迭代 算法 的 收敛 速 度 ， 获 得 高 精度 的 重 
构 结果 。 


~ 


图 2. Blob 重 构 模型 
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@ 背 投影 算法 (BPT) 


在 现 有 迭代 算法 中 ， 大 部 分 X2 的 选取 是 任意 的 ， 通 常情 况 下 选取 全 0 (或 全 1) 向 量 ， 
办 此 要 使 所 求 像素 值 X 达到 收敛 可 能 需要 迭代 多 次 。 为 了 能 够 减少 迭代 次 数 ， 加 快 算法 收 
敛 速度 ， 我 们 考虑 采用 直接 背 投影 算法 (BPT) 估计 像素 点 初 值 XIO。 直 接 背 投影 算法 是 
种 简单 的 重 构 算 法 , 它 认 为 每 个 像素 点 的 灰 度 值 是 所 有 通过 该 点 的 射线 投影 加 权 之 和 。 因 此 ， 
像素 点 灰 度 初 值 XO 可 表示 为 : 


M wp:; 
x = 下 
Wy 


然 背 投影 算法 估计 像素 灰 度 初 值 会 引入 星 状 伪 迹 , 但 比 任意 选取 初 值 的 方法 要 更 接近 
真实 的 灰 度 值 ， 从 而 能 够 加 快 收 剑 速度。 而且， 后 面 所 介绍 的 几 种 技术 能 够 有 效 地 抑制 星 状 
伪 迹 ， 降 低 误 差 引入 。 


@ 多 级 获取 策略 MAS 
现 有 和 迭代 算法 采用 的 是 逐 方向 更 新 策略 ,依次 顺序 对 相 邻 投影 方向 的 数据 进行 更 新 。 因 


为 顺序 投影 方向 之 间 存 如 


E 很 大 的 相关 


生 ， 容 易 引 入 误差 。 在 安排 


投影 数据 的 访问 方式 时 应 降 


低 相 邻 投影 数 


据 之 间 的 相关 怕 


生 ,使 各 个 相 邻 访问 的 投影 方向 尽 可 


能 地 正 交 ， 从 而 加 速 收敛 过 


程 。 为 了 满足 这 一 要 求 ， 我 们 采 月 


多 级 访问 策略 (MAS) 


重新 对 投影 方向 进行 排序 ， 使 一 


Pt 


系列 被 连续 访问 的 投影 角度 在 视角 范围 内 尽 可 能 地 均匀 分 布 。 己 经 证 明 当 两 个 方向 夹 角 为 
90°? 时 它们 的 相关 性 最 小 ， 而 第 三 个 方向 平分 前 两 个 方向 夹 角 时 三 个 方向 相关 性 最 小 。 基 于 
a ne 彤 方 癌 进行 章 新 排序 。 特 别 是 当 投 影 方 向 的 个 数 是 2 的 
窜 时 ,投影 方向 的 顺序 类 似 于 一 维 快速 傅 里 叶 变 换 。 假 设 投影 方 向 个 数 为 P， 每 个 方向 依次 
编号 为 0，1，……… ，P 了 -1。MAS 将 P 了 个 投影 方向 分 成 工 级, 工 可 以 表示 为 : L=[log,P]+1。 


在 第 一 级 ， 有 两 个 投影 方向 分 别 为 0” 和 90" ;， 接 下 来 的 第 二 级 也 包括 两 个 投影 方向 ， 
分 别 为 45” 和 135”; 第 三 级 有 四 个 投影 方向 ， 依次 为 22.5。、112.5。、675。 和 157.5。 
在 MAS 策略 中 ， 每 一 级 的 投影 角度 依次 平分 前 面 所 有 的 投影 角度 。 在 电子 断层 扫描 中 ， 
于 电镜 本 身 的 限制 ， 投 影 方向 无 法 覆盖 180”。 我 们 将 MAS 策略 进行 改进 ， 按 照 每 个 方向 
的 编号 进行 重新 排序 。 在 第 一 级 ， 选择 编号 为 0 和 2 的 投 晤 多方 向; 第 二 级 包括 编号 为 P/4 
和 (3/4) 书 的 两 个 投影 方向 ， 第 三 级 包括 编号 依次 为 P8、(5/8) P、(3/8) 了 和 (7/8) 了 的 
四 个 投影 方向 ， 后 面 的 各 级 以 此 类 推 。 在 最 后 一 级 时 ， 按 照 上 述 原则 计算 所 得 的 方向 可 能 
经 被 选择 了 , 我 们 需要 从 该 编号 位 置 开始 向 两 边 找到 距离 最 近 且 未 被 挑选 的 方向 。 在 迭代 
台 之 前 ， 所 有 的 投影 角度 按照 MAS 策略 进行 重 排 ， 然 后 每 次 迭代 时 按照 新 的 投影 角度 顺序 
进行 更 新 。 这 种 策略 能 够 大 大 降低 相 邻 方向 之 间 的 相关 性 ， 从 而 提高 算法 的 收敛 速度 。 


@ 自 适 应 误差 调整 ACC 
在 ASART 中 ,我 们 采用 一 种 数据 驱动 的 方式 对 x 进行 调整 : 第 磊 次 迭代 投影 误差 反 投 


影 补偿 时 不 仅 考虑 加 权 因 子 ， 而 且 还 考虑 x 中 自身 的 值 ， 所 以 第 i 条 射线 投影 对 像素 点 灰 度 
值 x 的 误差 补偿 为 : 


x ln -> Pi 


”Wo 


一 


h=1 


@ 列 集 合 替 代 技 术 CSS 
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在 ASART 中 , 我 们 采用 列 集合 奉 代 技术 CSS, 考虑 使 用 一 个 标量 8 表示 加 权 和 矩 阵 画 每 
列 和 的 最 大 值 的 倒数 。 而 矩阵 丈 中 每 列 和 的 最 大 值 可 以 认为 是 通过 一 个 像素 点 射线 的 最 大 
个 数 。 而 且 根据 丈 的 计算 模型 ， 通 过 一 个 像素 点 的 射线 的 最 大 个 数 即 为 投影 方向 的 个 数 。 

此 ， 此 标量 可 以 表示 为 : B=WB， 其 中 ，B 表示 投影 方向 的 个 数 。 采 用 CSS 技术 ， 达 代 


过 程 不 必 再 对 球 每 一 列 求 和 ， 可 以 大 大 节约 算法 运行 时 间 ， 同 时 也 能 节省 存储 空间 。 
根据 以 上 介绍 的 四 种 关键 技术 ， 改 进 后 的 ASART 算法 迭代 公式 如 下 所 示 : 


Dy wyP: ,j=1,2,...N 


-1 
Wy 
x = + p> 


wx Cp - ->， iDCN 


h=1 


(k) 


”om 


我 们 已 经 分 别 在 曙光 4000H、single GPU 和 multi-GPU 多 种 人 硬件 平台 上 实现 了 ASART 


算法 的 并 行 化 79。 


3.3 基于 GPU 平台 的 电子 断层 重 构 软 件 ATOM 


针对 现 有 电子 断层 成 像 技术 中 
重 构 软 件 的 不 足 , 特别 是 迭代 重 构 算 
法 速度 慢 的 缺点 , 我 们 开发 了 一 款 基 
于 图 形 处 理 器 (Graphic Processor 
Unit，GPU) 平台 的 电子 断层 重 构 软 : 


对 位 、 确 定 重 构 参 数 、 三 维 重 构 及 二 


于 : 在 二 维 对 位 方面 采用 了 迭代 的 平 
移 旋 转 对 位 方法 ， 提 高 对 位 精度 ; 在 
三 维 重 构 方面 ， 实 现 了 GPU 平台 上 : 
的 并 行 迭 代 重 构 算法 ， 其 中 Simulta- 
neous Iterative Reconstruction Tech- 
niques (SIRT) 算 法 达到 47 倍加 速 比 。 
为 保证 与 其 他 电镜 软件 的 融合 性 , 我 
们 对 软件 操作 过 程 中 的 电镜 序列 图 


获得 电镜 
二 维 图 像 序列 


件 ”ATOMP3]， 实现 了 二 维 投 最 交 图 像 es 


快速 建立 
重 构 模 型 


维 数 据 的 可 视 化 等 功能 。 其 特点 在 | -国信 站 移 


像 采 用 主流 的 MRC 格式 存储 。 图 3. ATOM 软件 基本 流程 
ATOM 软件 的 设计 流程 如 图 3 所 示 ， 其 中 实 线 框 内 为 ATOM 1.0 已 实现 的 功能 


AIOM 1.0 是 基于 Qt 和 CUDA 库 开 发 的 跨 
平台 软件 ， 使 用 C++ 和 C 语言 开发 完成 ， 其 中 
界面 代码 由 C++ 完成 ,核心 代码 主要 由 C 完成 ， 
通过 编译 可 以 运行 在 支持 Qt 和 CUDA 的 所 有 操 
作 系 统 上 。ATOM 1.0 设计 时 将 界面 和 可 执行 程 
序 分 离 ， 用 户 可 启动 ATOM 图 形 界面 操作 ， 也 
可 直接 执行 可 执行 文件 , 在 命令 行 下 能 完成 所 需 
功能 。 考 虑 到 软件 的 易 用 和 可 理解 性 ， 我 们 将 
AIOM 1.0 的 图 形 界面 设计 为 “基本 功能 ”“ 参 


49 


图 4 ATOM 界面 


数 交 
以 根 


http: 


心 
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互 ”、“ 运 行 信息 反馈 ”以 及 “显示 ”4 个 部 分 如 图 4)。ATOM 为 开源 软件 ， 用 户 可 
据 自 己 的 实际 需求 ， 添 加 相应 的 代码 ， 扩 充 相应 功能 。ATOM 的 下 载 地 址 : 


//feilab.ibp.ac.cn/software/atom/atom.html. 


目前 ATOM 已 完成 2.0 版 本 的 更 新 ， 图 中 的 虚线 框 部 分 的 功能 也 已 实现 。 


总 结 与 未 来 工作 


本 文 详细 介绍 了 电子 断层 三 维 重 构 技术 的 发 展 和 现状 , 分 析 了 当前 存在 的 主要 问题 , 着 
绍 了 我 们 在 电子 断层 三 维 重 构 方 面 的 相关 研究 工作 : 提出 一 种 新 的 基于 特征 提取 的 二 维 


Ce on a 
导 到 更 为 精确 的 重 构 结 果 ; 开发 了 国内 首 款 基 于 GPU 平台 的 电子 断层 重 构 软件 AIOM。 


在 下 一 步 研究 工作 中 , 我 们 将 主要 在 以 下 四 个 方面 继续 深入 开展 研究 : (1). 研究 新 的 无 


交 体 金 二 维 图 像 对 位 算法 ， 进 一 步 提 高 对 位 参数 的 精度 ; (2). 在 ASART 基础 上 ， 继 续 进 一 


步 提 


高 重 构 算法 的 收敛 速度 和 精度 ; (3). 开发 AITOM2.0， 研 究 在 GPU 集群 上 的 电子 断层 三 


维 重 


构 并 行 ，(4). 同 生物 学 家 合作 ， 将 电子 断层 三 维 重 构 技术 与 其 他 生物 成 像 技术 〈 如 光学 


PLAM 成 像 ) 相 结合 ， 开 展 光 电 联 合 三 维 成 像 的 研究 。 
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